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Abstract In these pedagogical lectures, we present the techniques of the AdS/CFT 
correspondence which can be applied to the study of real time dynamics of a 
strongly coupled plasma system. These methods are based on solving gravitational 
Einstein's equations on the string/gravity side of the AdS/CFT correspondence. We 
illustrate these techniques with applications to the boost-invariant expansion of a 
plasma system. We emphasize the common underlying AdS/CFT description both 
in the large proper time regime where hydrodynamic dynamics dominates, and in 
the small proper time regime where the dynamics is far from equilibrium. These 
AdS/CFT methods provide a fascinating arena interrelating General Relativity phe- 
nomenae with strongly coupled gauge theory physics. 



1 Introduction 



The current experimental program of heavy-ion collisions at RHIC and the forth- 
coming experiments at LHC open an interesting window onto properties of QCD 
matter at high temperatures, where it appears in the guise of a new phase — the 
quark-gluon plasma. At asymptotically high temperatures it should be a free gas 



Romuald A. Janik 

Institute of Physics, Jagiellonian University 
ul. Reymonta 4, 30-059 Ki'akow, Poland 
e-mail: romualdSth . if . u j . edu . pi 



1 



2 



Romuald A. Janik 



of quarks and gluons, however, at the experimentally accessible energies there are 
strong indications that the quark-gluon plasma is indeed a strongly coupled system 
(see e.g. 

This poses numerous problems for its theoretical description, yet at the same time 
makes its study theoretically interesting. One can roughly differentiate the physical 
properties of the quark-gluon plasma system into two broad categories — static and 
dynamic (real-time) properties. 

The first of these, the static properties, typically involve the study of equilibrium 
thermodynamics, the entropy, energy density as a function of temperature and more 
generally properties which can be directly deduced from the Euclidean formulation 
of finite temperature gauge theory. In this case lattice QCD is an effective tool for 
accessing these properties in the nonperturbative, strongly coupled regime. It deals 
directly with QCD and yields quantitative results directly applicable for the QCD 
quark-gluon plasma. 

The second class, the real time dynamic properties of strongly coupled plasma 
are much more difficult to access. They have to be formulated directly in Minkowski 
space and since lattice QCD methods are inherently Euclidean, it is very difficult to 
extrapolate numerical results to Minkowski signature. Moreover, it is exactly these 
kind of properties which are particularly relevant for the quark-gluon plasma pro- 
duced in heavy-ion collisions. 

To this end let us recall schematically the basic stages of a heavy-ion collision. 
First the two ultrarelativistic nuclei collide and the plasma is produced in a state 
very far from equilibrium. Then in a relatively very short time, it becomes ther- 
malized (or at least the pressure becomes isotropic with the residual anisotropy 
wholly due to flow). From that point on, hydrodynamic phenomenological models 
seem to describe the properties of the expanding plasma quite well |[2l- In partic- 
ular the plasma expands and cools, and when the temperature reaches the confine- 
ment/deconfinement phase transition one expects hadronization to occur 

It would be interesting to understand these various stages of the plasma dynam- 
ics directly from first principles. E.g. one would like to derive the hydrodynamic 
behaviour from some theoretical framework and not only use it as a phenomenolog- 
ical model. But what is even more interesting is the understanding of the process 
of thermalization and what governs the short thermalization time necessary for the 
applicability of hydrodynamic models. 

Unfortunately, for the case of QCD we lack appropriate theoretical methods 
which would be applicable to these kinds of problems at strong coupling. A possible 
route that one may take is to consider an analogous set of problems in a different 
theory for which appropriate real-time nonperturbative tools exist. 

The new method for studying nonperturbatively various gauge theories (although 
not directly QCD) is the AdS/CFT correspondence |3|, which translates dynami- 
cal problems in strongly coupled gauge theory into (usually) gravitational ones in 
higher number of dimensions. Its main advantage is that it works equally well in 
Minkowski as well as in Euclidean signature. In these lectures we will consider the 
AdS/CFT correspondence in its simplest original setting for the maximally super- 
symmetric — 4 Super- Yang-Mills theory. Of course, one has to be aware that the 
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cost of switching the theory of interest from QCD to jV — 4 SYM is that we may 
most probably loose direct quantitative applicability of our results to realistic heavy- 
ion collisions. Moreover, there are certain marked differences between the theories 
which will cause a huge difference for certain physical phenomenae, at the same 
time being unimportant for other questions. We will discuss some of these points in 
these lectures. 

Nevertheless, we would like to point out that currently we do not have any gauge 
theory in which we would have a theoretical understanding of the issues described 
earlier. Therefore it is very interesting to study these issues for the case of jY = 4 
SYM and use the results as a point of reference for analyzing the situation in QCD. 
Later, one could try to generalize these results to more complicated versions of the 
AdS/CFT correspondence for theories closer to QCD. In essence, this motivation is 
in line with the statement that jY ~ 4 SYM is the harmonic oscillator of four di- 
mensional gauge theories. If one tries to develop some theoretical tools, one should 
better first apply them to the 'harmonic oscillator'. 

We have tried to make these lectures very pedagogical and self-contained. Our 
main emphasis in the presentation is to show how one can use the AdS/CFT cor- 
respondence as a tool even for far from equilibrium configurations without presup- 
posing any kind of dynamics (which are in fact not known for nonlinear far from 
equilibrium systems). Therefore our presentation of hydrodynamics is subordinate 
to this goal, especially as a very general discussion focused on hydrodynamics per se 
is contained in the lectures by V. Hubeny at this school |4|. 

The plan of these lectures is as follows. In section[2] we introduce the AdS/CFT 
correspondence, in section [3] we compare some properties of plasma in the ,yY — 4 
theory and in QCD. Then we proceed to present the AdS/CFT setup specialized to 
the study of time-dependent dynamics of strongly coupled plasma. We then illustrate 
these methods in section [5]by discussing two important examples - the appearance 
from this setup of the standard planar AdS black hole, and a planar shock wave. 
In this section we also discuss some subtleties arising with different choices of co- 
ordinate systems which will be relevant later. Then in section |6j we introduce the 
main physical example of a time-dependent plasma configuration — the boost in- 
variant flow. In section [7] we analyze its large proper time asymptotics and show 
how nonlinear perfect fluid dynamics arises from the AdS/CFT methods. In section 
[8] we show how one can see first corrections coming from shear viscosity, and in 
section|9] for completeness, we will summarize briefly the current status of hydro- 
dynamics in AdS/CFT. Then in section 10 we introduce physical situations, where 
the plasma dynamics is not describable by hydrodynamics, and finally, in section[TT| 
we apply the AdS/CFT methods to study boost invariant flow in the far from equi- 
librium small proper time regime. We close the lectures with conclusions and an 
appendix with a short guide to the literature regarding work done on related topics 
which were not covered here. 
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2 The AdS/CFT correspondence 

The AdS/CFT correspondence fi\, in its original form states the equivalence of two 
apparently completely different theories: the yy = 4 supersymmetric Yang-Mills 
theory (SYM) in four dimensions and type JIB superstrings in a ten-dimensional 
curved AdSs x background. Since then, it has been generalized in various di- 
rections, extending it to a wider class of gauge theories, at the cost of making the 
dual string backgrounds more complicated. Throughout these lectures we will stay 
within the context of the AdS/CFT correspondence for ^ = 4 SYM theory. 

The reason why the AdS/CFT correspondence is so interesting is that the non- 
perturbative strong coupling regime of the ,yK = 4 gauge theory is mapped to the 
(semi-)classical strings or just (super)gravity which, in contrast to the gauge theory 
side, is at least theoretically tractable. Therefore one can use the AdS/CFT corre- 
spondence as a new method for accessing the very difficult nonperturbative gauge 
theory physics. 

The AdS/CFT correspondence is an equivalence, so in principle any state/pheno- 
menon on the gauge theory side should have its direct counterpart on the string side 
and vice-versa. However one should keep in mind that the correspondence is an 
equivalence of gauge and string theory, so the dual counterpart does not have to be 
in the well understood (super)gravity sector. Fortunately, it will turn out that for the 
questions considered in these lectures namely the study of the dynamics of plasma 
expansion, the dual description will be purely gravitational. 

Apart from the direct physical interest, the AdS/CFT correspondence is also the- 
oretically very interesting as it translates various dynamical gauge theory questions 
into a geometrical language described by higher-dimensional General Relativity 
(GR). This leads to quite fascinating links between the two fields, providing a whole 
range of physically motivated interesting questions which could be addressed by GR 
methods. In the other direction, various notions introduced by the GR community 
like dynamical apparent horizons find their application and new interpretation on 
the gauge theory side. 



Effective degrees of freedom at strong coupling 

As an illustration of the use of the AdS/CFT correspondence, and as a justification 
for the gravitational methods let us consider the question of finding effective degrees 
of freedom for strongly coupled ^ — 4 SYM. By the AdS/CFT equivalence, it 
amounts to asking the same question for superstrings in AdS 5 x S^. 

Let us first recall the case of closed strings in flat space. The string worldsheet 
action is characterized by a dimensionfull parameter a' (related to the string ten- 
sion). The various vibrational modes of the string correspond to particles (fields) 
with distinct masses ^ 

- ^ (1) 
a' 
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The massless modes correspond to the graviton (and its whole supergravity mul- 
tiplet). There is also an infinite tower of massive modes. In the limit of a' 0, 
the massive modes become very heavy and effectively decouple at fixed energies 
leaving as the governing dynamics just (super)gravity. 

In the case of strings in AdSs x , the a' parameter becomes proportional to 
1 / \fX where A = gyM^c '^^e 't Hooft coupling of the dual = 4 gauge theory. 
The vibrational modes again split into a massless (super)graviton multiplet and a set 
of massive modes. The formula ([T]i is no longer exact, but the parametric behaviour 
with a' still holds. So in the strong coupling limit, those massive string modes be- 
come very heavy and effectively decouple leaving essentially supergravity modes as 
the effective degrees of freedom. Since the AdS/CFT correspondence postulates an 
equivalence with gauge theory, these should also correspond to the effective degrees 
of freedom of the gauge theory at strong coupling. 

Once we lower the coupling, the massive modes become lighter and their effects 
will no longer be negligible. Initially, their effects may be absorbed into corrections 
to the gravitational Einstein-Hilbert action (so-called a' corrections), but at low 
coupling corresponding to the perturbative regime the spacetime description is not 
known. 

Finally, let us note that the ,yV — 4 SYM theory is quite special in that it allows 
for such a clean separation between gravity modes and massive string modes. Pre- 
sumably a dual description of real QCD (or even of large A^^- pure YM) would not 
have such a property. 



3 Why study ,yV = A plasma? 

Since we will be using the AdS/CFT correspondence for ^ — 4 SYM as a calcu- 
lational tool for analyzing strongly coupled dynamics of gauge theory plasma, we 
will be essentially considering plasma in the supersymmetric o/K ~ 4 gauge theory. 
This theory is completely different from QCD at zero temperature. It is supersym- 
metric, exactly conformal, does not have confinement. However once we turn on 
some nonzero temperature (or consider not the vacuum but some appropriate state), 
supersymmetry is broken and temperature (or energy density) introduces a scale. 
So qualitatively, we may expect to have similarities with QCD plasma in a certain 
window of temperatures where it is strongly coupled, approximately conformal and 
(by definition) deconfined. 

However we have to keep in mind some definite differences w.rt. QCD plasma. 
Firstly, the ,yV = 4 theory has no running coupling so, in contrast to QCD, even at 
very high tempeartures/energy densities the coupling may remain large. Secondly, 
the equation of state of the o/K = 4 plasma is exactly conformal {E = 3p) which is 
only an approximation for a certain range of temperatures for QCD plasma. Indeed 
we know, from lattice QCD, that deviations from a conformal equation of state 
are important close to 7^. In addition, we have other consequnces of conformality 
like that the bulk viscosity for the c/K ~ 4 theory is exactly zero. Thirdly, for the 
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jV — 4 theory (on Minkowski spacetime) there is no phase transition — no analog 
of the confinement/deconfinement phase transition of QCD. Therefore as the plasma 
expands and cools, in the jV ~ 4 theory it will expand indefinitely, while in QCD it 
will cool down to the phase transition temparature and hadronize. 

From the above discussion we see that the applicability of using jV = 4 plasma 
to model real world phenomena depends on the questions asked. It may give a good 
qualitative picture for the range of temperatures where we have similarities. How- 
ever, let us note that in this theory we may compute the dynamics from 'first princi- 
ples' (using the AdS/CFT correspondence). For QCD, unfortunately, we do not have 
any similar calculational technique, even numerical, which would enable us to study 
real-time dynamics of the strongly coupled quark-gluon plasma. Therefore it is in- 
teresting to build up results on strong coupling properties of ~ 4 plasma and use 
them as a point of departure for analyzing or describing QCD plasma. Eventually, 
one might consider more realistic theories with AdS/CFT duals which are closer to 
QCD. In those cases generically the dual gravitational backgrounds are much more 
complicated so it is advantageous to start from the simplest setting for the jV = 4 
SYM theory. 

Another motivation for studying the dynamics of jV ~ 4 plasma is that the natu- 
ral language of the AdS/CFT correspondence is quite new w.r.t. conventional gauge 
theory methods. So by studying relatively simple examples we may build up some 
new physical intuitions within this novel language. Also the interrelations with Gen- 
eral Relativity physics are fascinating from the purely theoretical point of view. Last 
but not least, there may be some unexpected discoveries like the celebrated univer- 
sality of the shear viscosity to entropy ratio r\ls |5|, which, at strong coupling, 
remains equal to 1 /4;r for any theory with a gravitational dual L6J (see the lectures 
by A. Starinets at this school Q). 



4 The AdS/CFT setup for studying real-time dynamics of plasma 

Let us now briefly review the AdS/CFT correspondence on a more technical level, 
concentrating on the features relevant to the study of the time evolution of a plasma 
system. 

The factor in the AdSe, x 5^ background is associated with a global 5(9(6) = 
5t/(4) symmetry of the jy = 4 theory. In the following, we will only consider 
systems which do not break this symmetry, so the 5^ factor will be irrelevant and 
the whole dynamics will be concentrated in the AdS^ factor 

The 5-dimensional Anti-de-Sitter spacetime AdS^ can be given by the following 
metric 

2 rjavdxf^dx^ +dz^ 
ds = 2 (2) 

with z > 0. z = is the boundary of AdS^, while the region z > is often called 'the 
bulk' . This choice of coordinates covers the Poincare patch of global AdS^ and is 
relevant for the case when the dual gauge theory lives in M''^ Minkowski spacetime. 
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The above geometry can be understood to correspond to the gauge theory vacuum 
state. In particular gauge theory operators Uke the energy-momentum tensor T^y all 
have vanishing expectation values in this state 

(7mv)=0 (3) 

Let us recall that the AdS/CFT correspondence states the equivalence with super- 
strings in AdSsi xS^). So, on the string side, we may excite any normalizable mode, 
in particular we may excite gravitons. This will correspond to some states in yy = 4 
SYM with (7}iv) 7^ 0. We expect a configuration of — 4 plasma to be a very 
complicated state which would correspond to exciting very many gravitons. Then it 
is better to interpret it instead as a change of the background: 



ds^ = gl''pdx"dxP = ' ' ^2 (4) 

with the metric coefficients being now generic functions of all the five coordinates. 
We therefore seek to describe a plasma configuration in terms of the geometry 
gtivix^ ,z). Let us note that the above choice of the metric Q is always possible 
after a suitable change of coordinates. Such coordinates, in which the metric has the 
form (|4]) are called Fefferman-Graham coordinates. 

The geometry Q cannot be, however, completely arbitrary. It must form a con- 
sistent background for strings, so it must satisfy 5-dimensional Einstein's equations 
with a negative cosmological constanj^ 

Rap-lgl1iR-6gfp=0 (5) 

Furthermore, for a physical state in the gauge theory this geometry should not have 
a naked singularity. This turns out to be a crucial requirement with far reaching 
consequences for the resulting dynamics of the yy — 4 plasma, as we will see in the 
following. 



The gravity — > (T^v) dictionary 

Once one has the geometry Q corresponding to some plasma configuration, the 
key question is what is the energy momentum tensor of that gauge theory system. 
The way to derive the answer, holographic renormalization, has been explained in 
the lectures by K. Skenderis [Sj. Here we just summarize the outcome derived in 
||9l, which in the Fefferman-Graham coordinates defined by Q takes a particularly 
simple form. 



These equations are equivalent to the original ten-dimensional type IIB supergravity equations 
when we preserve full S0(6) symmetry of and no other fields are turned on. The negative 
cosmological constant is a remnant of the RR 5-form under this dimensional reduction. 
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Suppose that the metric coefficients g^v{x^ ,z) have the following Taylor expan- 
sion near the boundar\n 



The spacetime dependence of the energy-momentum tensor carries a lot of informa- 
tion about the dynamics of the plasma - its energy density, momentum flow, stress 
tensor Indeed, it is just this information which is exactly the direct outcome of hy- 
drodynamic simulations of realistic heavy-ion collisions. Finally, let us emphasize, 
to avoid any chance of confusion, that the T^y is the energy-momentum tensor of the 
dual gauge theory. On the gravity side we are always dealing with vacuum Einstein's 
equations. 

The construction outlined above leads to the following scenario of investigat- 
ing a plasma system in strongly coupled ^ — 4 SYM. One starts from some ini- 
tial conditions for the 5-dimensional Einstein's equations. Then the geometry is 
evolved forward in time by solving Einstein's equations. Finally using the above 
formula (jTji, one extracts the (7/iv) of the corresponding plasma system. The details 
of the evolution of T^y are very interesting from the point of view of physics, espe- 
cially in the far from equilibrium region, where very little is known about thermal- 
ization/isotropisation and transtion to a hydrodynamic expansion. We will follow 
this route in section [TT] 

The (Tfjiv) — > gravity dictionary 

It turns out to be very fruitful to ask also the opposite question. Suppose that we are 
given a certain spacetime profile of the energy momentum tensor (7)jv) - how to 
construct the dual 5-dimensional geometry? The prescription is really just running 
the preceding recipe backwards. One has to solve Einstein's equations 



(6) 



Then the expectation value of the energy-momentum tensor is 




(7) 



1 50 n 



= 



(8) 



with the boundary condition 



(9) 



where glilix^) is related to (^T^vix^)) through 



^ Here we always assume that the gauge theory hves in fiat Minkowski space, hence the leading 
J]^v and the absence of a term (see (5) for details). 
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It turns out that for a solution to exist, g\iv{xP) has to be traceless and conserved, 
which is of course physically expected for an energy-momentum tensor in a confor- 
mal theory. However here it is just an independent consequence of 5-dimensional 
Einstein's equations. 

Once such a ^|fy {xP ) is chosen, Einstein's equations determine uniquely the so- 
lution in the bulk (at least locally i.e. all higher coefficients of the Taylor expansion 
of gfiv {x^ , z) around z = are uniquely determined). 

In this way we see that for every energy momentum profile which does not violate 
the standard requirements of energy-momentum conservation and tracelessness we 
may construct a dual gravity background. However generically, such a geometry 
will be highly singular with naked singularities in the bulk. The requirement of 
the absence of naked singularities will very strongly constrain the admissible bulk 
geometries and hence also the possible spacetime profiles of the energy momentum 
tensor This is a nontrivial constraint on the dynamics of the gauge theory as the 
spacetime profile includes of course the time evolution of 7j,v 

This line of reasoning was introduced in [lOJ as a way of determining the pos- 
sible evolution of the energy momentum tensor. One first picks a family of profiles 
Tjjv{xP), then one constructs for each of them a dual geometry by solving Einstein's 
equations with appropriate boundary conditions. Finally, one picks the allowed dy- 
namics by requiring that the corresponding dual geometry is nonsingular. We will 
describe this procedure in detail in the first part of these lectures. 



5 Exact analytical examples 

In this section we will illustrate the AdS/CFT methods by analyzing two simple 
examples of plasma configurations for which the dual gravitational background can 
be computed in closed form ifTOl . Both of these examples have also a clear physical 
interpretation. 



5.1 A case study: static uniform plasma 



Let us first consider the simplest configuration of plasma, namely with a uniform 
and static distribution of energy density filling up the whole spacetime. The energy 
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momentum tensor is just a constant diagonal one: 



Tjxv — 



/£ 0\ 
0/700 
0/90 



(11) 



with £ = 3/7. In order to find the dual gravity background, we have to solve Ein- 
stein's equations with the boundary conditions given by (j9])-([T0]). Due to the fact 
that the energy momentum tensor is constant, the metric will only depend on the 
z coordinate and the Einstein's equations reduce to ordinary differential equations 
which can be solved explicitly. The result is 



ds"- = - 



dt^ + {l+zVzl)-^ + 



,2 



where the parameter zq is related to E through 



E = 



3Nl_ 



(12) 



(13) 



Although it is not evident at first glance, the geometry (12 1 is exactly the standard 
AdS planar black hole 1 11 1, but written in the Fefferman-Graham system of coordi- 
nates. We will give the explicit form of the coordinate transformation to the standard 
AdS Schwarzschild form shortly. 

The fact that the dual geometry turns out to be a black hole has significant impli- 
cations for the physics. Let us note that we did not assume that a black hole would 
appear It came, in a unique way, from solving Einstein's equations with appropriate 
boundary conditions. 



The parameter zo appearing in (12 1 is the location of the black hole horizon. The 
Hawking temperature Th which is given by 



Th 



V2 



(14) 



is identified with the gauge theory temperature. This may be most easily seen by 
computing the Hawking temperature through a Wick rotation of the metric ( [T2| and 
requiring the absence of a conical singularity at z = zq. This requirement leads to a 
specific periodicity condition for the Euclidean time coordinate which is inversely 
proportional to the Hawking temperature. But, according to the AdS/CFT corre- 
spondence the metric induced on the boundary z = (up to an overall rescaling 
by z^) is exactly the metric of the (now Euclidean) gauge theory. Thus the gauge 
theory also has a compactified Euclidean time with the radius given by the same 
temperature ifTTl . 

Another gravitational concept which carries over to the dual gauge theory is the 
Bekenstein-Hawking entropy which is identified with the entropy of the dual gauge 
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theory plasma system. In this case, the entropy per spatial 3-volume is 

Now we can use the relation between the horizon parameter zq and temperature to 
express the result completely in terms of gauge theoretical quantities 

s = ]fh''T^ (16) 

Finally, we may use the relation between the energy density E and zo, and the 
link with temperature to express the energy density as a function of T . We get 

E = In^k^T'^ (17) 
8 

The nontrivial factor here is the numerical coefficient, which is different from the 
corresponding one for the free massless gas (Stefan-Boltzmann). Similarly, the en- 
tropy density derived earlier is 3/4 of the free gas answer lfT2l . This mismatch is 
quite natural since here we are dealing with a strongly coupled plasma. In fact 
similar deviations from the Stefan-Boltzmann answer have been observed in lat- 
tice studies of QCD thermodynamics above the confinement/deconfinement phase 
transition. 

Before we move on to discuss various systems of coordinates for this geometry, 
let us note that it is exactly this geometry which is used to describe ^ — 4 SYM 
at fixed nonzero temperature T. This interpretation is obvious from the above men- 
tioned Euliclidean continuation, but can also be understood directly in Minkowski 
signature, where a hnk with the real-time formalism of finite-temperature QFT ap- 
pears lfT3l[T4l[T5][T6l. 



Various coordinate systems 

The geometry ( [T2] i has been presented in the Fefferman-Graham coordinates, in 
which the connection to the gauge theory energy-momentum tensor is simplest. 
However these coordinates have also some significant drawbacks, of which one has 
to be aware. 



Let us first perform a coordinate transformation to bring the metric ( 12 1 into the 
standard AdS Schwarzschild form. To this end set 

Zs,d = I ^ (18) 

^i+^V4 
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Then, the metric becomes 



ds- = ^dt' + ^+ ^ — (19) 

^std ^std ^ ^stdl^Ql ^std 



with zo = zo/\/2. Looking at the transformation of coordinates ( 18 1, we see that the 
Fefferman-Graham coordinates cover only the region between the boundary and the 
horizon. Even if one would analytically continue the metric for z> zq one does not 
go beyond the horizon but rather returns back to the boundary. 

The standard Schwarzschild coordinates also break down at the horizon and, in 
order to have explicit regularity at the horizon, it is convenient to introduce yet 
another system of coordinates — the (ingoing) Eddington-Finkelstein coordinates. 

These coordinates may be obtained from the standard Schwarzschild ones by 
redefining the time coordinate: 

. 1 - ^ Zstd , , ZO + Zstd \ 

tEF = t- -ZO 2 arctan h log (20) 

4 V 20 ZQ-ZstdJ 

The metric becomes then 

1 ^^zt^/zn 1 dtppdzad dx} 
ds^ = f^dtlp+2 "'' +-^ (21) 

^std ^std ^std 

The crucial advantage of these coordinates is that the horizon is a perfectly regular 
point and one can enter the region inside the horizon. The lines = consff^ are null 
geodesies falling into the black hole and reaching the singularity at Zstd = °°- These 
coordinates were used extensively in V. Hubeny's lectures at this school with Zstd 
substituted by r = I / Zstd, which brings them to the canonical form. 

Finally let us note that in the formula ( [20| , the time coordinate gets an infinite 
shift at the horizon. In the case of the static black hole geometry this is completely 
harmless as the metric is time-independent, however for the time dependent geome- 
tries which will be the focus of these lectures this shift will give rise to some spuri- 
ous singularities in the Fefferman-Graham treatment (fortunately appearing only at 
NNNLO in the large proper time expansion of the geometry). 

Some comments are in order here. Of course in General Relativity nothing de- 
pends on the choice of coordinate system. This is true if we are dealing with an 
exact solution of Einstein's equations - we may analyze it in any coordinate system 
we like. However if we perform an expansion in time in some coordinate system 
and deal with approximate solutions truncated at some order, we may get spurious 
singularities like in Fefferman-Graham at third order 

In these lectures we will nevertheless present the analysis in Fefferman-Graham 
coordinates (apart from a brief review of the Eddington-Finkelstein approach of 
ifTTl in section |9]l. The general formulation in Eddington-Finkelstein is considered 
in detail in the lectures of V. Hubeny, and has been applied to the boost invariant 
setting in 1191 and EOl |2TI . One motivation for this choice of presentation is that 
our main focus is in reaching the small proper time regime, where we deal with 
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exact solutions of the Einstein's equations and hence do not need to worry about 
these subtleties. Also there, the analysis of initial conditions in Fefferman-Graham 
coordinates is simpler 



5.2 A case study: a planar shock wave 



Another case of a gauge theory energy-momentum tensor for which the dual geom- 
etry is exactly solvable is a planar Shockwave concentrated on the boundary. The 
only nonvanishing component of T^y is 

r__=M/(x-) (22) 

Such a configuration, for f{x^) — lJ.5{x^) represents a planar shock wave of gauge 
theory matter moving at the speed of light along one light cone direction. It may be 
understood to represent an analog of an ultrarelativistic nuclei. Then the dual metric 
is found to be 

T dx'^dx^ + f(x^)z^dx^ +dx^,+dz^ 
ds^ = — — '—^ ± (23) 

This configuration was proposed in fTol, being the simplest member of a family of 
shock wave solutions with.JCx dependence derived in |22| (see also l23l ). A natural 
question to consider is a collision of two such shock waves, one propagating along 
the light cone direction, the other along x'^ . There have been some preliminary 
investigation along these lines in 1 24 25 , 26 1 . However, a complete analysis remains 
to be done. Finally, let us note that this kind of shock wave is sourceless in the bulk, 
in contrast to the shock waves considered recently in ll27ll and in the lectures of 
A. Yarom ||28| . For references on work done on these other kinds of shock waves 
consult im. 



6 Boost-invariant flow 

Let us now concentrate on a concrete evolving plasma system. Since eventually we 
would like to solve exactly Einstein's equations, one has to introduce as much sym- 
metries as possible to reduce the complexity of the task, at the same time allowing 
for nontrivial physics to intervene. A natural choice in the context of heavy ion 
collisions is the requirement of longitudinal boost invariance. This assumption was 
introduced by Bjorken 1291 back in 1983 to model ultrarelativistic collisions. Basi- 
cally, the motivation is that at infinite energy, a finite boost along the collision axis 
would not modify the physics. This is certainly not an ideal approximation, however 
it is used in basically all hydrodynamic codes for modelling relativistic heavy-ion 
collisions at RHIC I2j|. We will make here an additional assumption that there is 
no dependence on the transverse coordinates, which corresponds to the limit of in- 
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finitely large nuclei. This is not really necessary for discussing the hydrodynamic 
limit (see HTHH) but will be essential in the far from equilibrium regime of small 
proper times. 

When assuming boost invariance, it is natural to pass to proper- time/spacetime 
rapidity coordinates {z,y,xi.X2) 

t = Tcoshy JC3 = Tsinhy (24) 

Then it turns out that the only non-vanishing components of the energy-momentum 
tensor are T^^, r„ and r„ = T^^xi = Txxn- Moreover, these components become 
functions of T alone. 

We should now impose tracelessness =0 and conservation of energy momen- 
tum T-y^ = condtions, which take the form 

-Trr + \Tyy + 2T„ = 
d 1 

dx ■' 
These equations determine 7|iv uniquely in terms of a single function e(T) 



/e(T) \ 



-x^-^e{T:)-x^e{T) 

1 ^ d 



£(T)+iT^e(T) 

V e(T)+iT^e(T); 



(25) 



The remaining function e(T) can be interpreted as the energy density of the plasma 
at mid-rapidity (i.e. at xt, = 0) as a function of (proper-) time. 

Let us note that the above decomposition was purely 'kinematical' — valid in any 
conformal 4D theory at any coupling. The determination of £(t) will be an issue 
of understanding the dynamics of the theory of interest — here = 4 SYM. In 
particular, weak coupling perturbative considerations I.30J lead to the free streaming 
behaviour 

e(T) ~ ^ (26) 

If, on the other hand, we would suppose that the plasma system behaves as a perfect 
fluid, then on top of the decomposition ( [25] l we would impose 

T^''' ^{e + p)u^'u'' -pn^" (27) 



with e = 3p. By our symmetry assumptions m'' = (1 , 0, 0, 0) and we get in particular 

' Tyy 



p = \Tyy — r„, which gives a differential equation for e(T) 



T^e(T)-e(T) = e(T)+iT^e(T) (28) 
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with the celebrated Bjorken solution 

, , const. ^^^^ 
e(T) = (29) 

Other dynamical assumptions would modify the functional form of e(T). E.g. if the 
fluid would not be a perfect fluid but would have a nonzero viscosity (proportional 



to r ~' e 4 as should be the case for a conformal theory), then ( 29 1 would no longer 
be exact but would have corrections starting with 

+ (30) 

with J7o related to the shear viscosity. Here we set a single dimensional scale to 
unity. It can be easily reinstated in all terms by dimensional analysis. Further 1 / t^/^ 
corrections are uniquely determined in terms of rjQ. If the dynamics would follow 
2nd Qj.(jg]- viscous hydrodynamics, these 1 /t'*/^ corrections would be different and 
would involve additional, second order transport coefficients. 

So we see that the knowledge of e(T) contains a lot of information on the dy- 
namics of plasma. In the rest of these lectures our goal will be to deduce what e(T) 
is singled out by the AdS/CFT correspondence. Initially we will concentrate on the 



large T asymptotics of e(T) and then, in section 1 1 move to consider the behaviour 
of e{z) for smaU t. 



7 Large proper time behaviour 

Let us first concentrate on the large T asymptotics of e(T) and consider determining 
the exponent s in 

£(t) — forr^oo (31) 

We win follow here the strategy outlined in section]?] and construct, for each s, the 
dual geometry. Then we will check for which s the dual geometry is nonsingular. 
This condition will determine s. This approach was proposed in [10], where more 
details may be found. But first, let us narrow down the range of possible s. We will 
demand that the energy density in any reference frame is non-negative i.e. 

Tfjivt'^t" > (32) 

for any timelike 4- vector t^. This leads to the inequalities 

£(t)>0 e'(T)<0 Te'(T) > -4£(t) (33) 

In particular, for ( [3T) , we obtain that < s <4. 
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Now we have to construct the dual geometry to a plasma configuration with the en- 
ergy density behaving like ( pT) . Since the geometry will have the same symmetries 
as assumed for the plasma system, we are led to the following ansatz 

^,2 ^ ^ [-e<^^')dT^+e''^^''hW + e'^^^''Uxl) + ^ (34) 

Again, let us reiterate that the above ansatz is completely general. At this stage the 
choice of the Fefferman-Graham system of coordinates is perfectly legitimate. 

According to the approach explained in section |4j we have to solve Einstein's 
equations 

^a/3-^ga>-6g'«^=0 (35) 
with the boundary conditions 

a(z,T) = -z^e{x)+z^a(,{x) +z^ai{x) + . . . (36) 

It is instructive to find the explicit form of the first few coefficients in the above 
Taylor serie^ Using a computer algebra system we obtain 

+ 16 ^ ^+128t3 128t2 64t 384^ ^^ij ^ + ^^') 
Let us now specialize to the case of interest e(T) — 1/t^ We get 



where the terms dominant for large T are outlined in bold. Looking at the above 
formula and analyzing a couple of higher order terms we may convince ourselves 
that the dominant terms in a„(T) for large T will be of the form 

z I z 



z"a„{T) = — for large T (38) 

This shows that it is natural to introduce a scahng variable 



' This is an exact result without any approximation. 
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v=A- (39) 
Consequently, the metric coefficients will have an expansion of the form 

a{z,z)=ao{v) + ^ai{v) + ... (40) 

Several comments are in order here. Firstly, the appearance of the scaling variable 
at late times is a dynamical consequence of the structure of Einstein's equations. 
We will find later, that for small proper times a similar structure will not appea^ 
Secondly, the separation of dynamics into a scaling variable and an expansion in in- 
verse powers of T corresponds to a gradient expansion (c.f. 1 17 1 and the lecture by V. 
Hubeny |4|). Finally, the appearance of a scaling variable reduces the very compli- 
cated Einstein's equations to a system of nonlinear ordinary differential equations: 

v(2fl'(v)c'(v) + a'{v)b'{v) + 2b'{v)c'{v)) - 6a' {v) - 6b' (v) - 12c' (v) + vc'{v)^ = 
3vc'{vf + vb'{vf + 2vb"{v)+4-vc"{v) - 6b' (v) - 12c' (v) + 2vb' {v)c' (v) = 
2vsb"{v) + 2sb'{v) + 8fl'(v) - vsa'{v)b'{v) - 8fo'(v) + vsb'{vf+ 

+Avsc"{v) +Asc' {v) ~2vsa' {v)c' {v) +2vsc' {vf = 

which can be solved exactly. The solution is 

fl(v) =A(v) -2m(v) 

b{v) =A(v) + (25-2)m(v) 

c(v) =A{v) + {2~s)m{v) 



where 



with 



A(v) J- (log(l +A{s)v')+\og{l-A{s)v')) 
'«(v) (log(l + 4 (.) v^) - log(l - 4 (.) v^)) 



/3s2-8i + 8 

A{s) 



24 

Now we can analyze the singularities of the above geometry. We see that there is 
a potential singularity where the argument of the logarithm vanishes. Of course, it 
may well be a coordinate singularity, so we have to evaluate a curvature invariant 
like 

^H2=/?A:va/3/?^,„^ (41) 

Moreover, since our geometry is only exact in the scaling limit, we have to evaluate 
in the same limit i.e. T — > 0°, z — > 0°, keeping the ratio v = ^ fixed. The resulting 

* For the case e(T) — > const as T — > 0. 
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expression is quite complicated and can be found in ifTOl . Its general structure is 

Name rat or(v;s) 



(l-4(i)2v8)'* 



(42) 



We see that for generic s, there is a 4'^' order pole singularity in the curvature. It 
turns out that this singularity gets cancelled by the numerator only for a single value 
of s: 

4 

.= 3 (43) 

which is just the asymptotic scaling characteristic of perfect fluid hydrodynamics. 
In this way we see that nonlinear perfect fluid hydrodynamics arises at late stages 
of plasma expansion as a consequence of the AdS/CFT correspondence. 




hoiizon 



Fig. 1 The difference between using Fefferman-Graham and Eddington-Finkelstein coordinates 
for checking nonsingularity represented by arrows. 



Note that in this way we do not approach the 'horizon' directly, but approach it 
asymptotically along a 'parallel' trajectory. An analogous analysis using Eddington- 
Finkelstein coordinates enables us to pass through the 'horizon' and require directly 
the nonsingularity of the metric there. The difference between the two procedures 
is summarized in Fig. 1. Both methods give equivalent results up to two sublead- 
ing orders in the large proper-time expansion. In order to go beyond that, however, 
one has to use Eddington-Finkelstein coordinates (currently the only result in this 
direction beyond 2'"' order is in [31] in the boost-invariant setting). 



7.2 Perfect fluid geometry 

Let us now examine more closely the dual geometry corresponding to the perfect 
fluid value of s = 4/3. Then the complicated formulas for the metric coefficients 
involving generically irrational powers and square roots obtained above simplify 
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drastically and we obtairj^ 



£0 

3 t4/3 



3 t4/3 



{x^dy^ + dx\ 



dz^ 



(44) 



This geometry is analogous to a black hole (c.f. ([12])) with the position of the horizon 
moving into the bulk as 

zo = \7 - ■ (45) 
V eo 

This has a clear physical interpretation. Recall that for a static black hole, the po- 
sition of the horizon in the bulk is inversely proportional to the temperature. Thus 
here we have a dual counterpart of the plasma undergoing cooling during expan- 
sior0 Indeed naively generalizing the static formulas leads to 



T = 



-T 3 



^20 7Z34 



(46) 



A more indepth analysis of these phenomenae using the framework of event f32\ or 
dynamical |20 , 33 , 31 1 horizons has been made, although a complete understanding 
of the notions of temperature and entropy in the fully dynamical case is still lacking. 



8 Plasma dynamics beyond perfect fluid 

One of the key discoveries of the AdS/CFT correspondence was the derivation of 
the universal value of shear viscosity for plasma at finite nonzero temperature. This 
was done using linear response theory in |5|. It is thus interesting to examine how 
do viscous effects manifest themselves in the current nonlinear setting. 

Let us first examine the question whether we can see if the perfect fluid dynamics 
is violated from the dual gravitational point of view. Suppose that it is not and that 
consequently 

e(T) - \ (47) 

is valid at all proper times without any corrections. Then one can find the next orders 
in the large T scaling expansion of the metric coefficients 

a(z,T)=ao(v) + 4r«2(v) + ... (48) 

T3 



^ We reinstated liere a trivial dimensional scale eo- 

* The dual counterpart of cooling was first suggested qualitatively in 1341 . 
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This can be done explicitly since, fortunately, the equations for the corrections are 
linear albeit still quite complicated. Peforming this calculation, and computing the 
curvature up to this order yields 

m^=RapysR"^^^ = W +^^+... (49) 
nonsingular singular! 

where the indicated singularity is of the very strong 4''' order pole type. This means 
that the perfect fluid behaviour (47 1 has to be violated. 

Let us now be completely agnostic about viscous hydrodynamics and assume a 
completely generic type of corrections: 

eW^^fl-^l (50) 



T3 



Solving the Einstein's equations with the appropriate boundary conditions set by pOf , 
computing the curvaturaj yields 

nonsingular nonsingular singular] singular! 

The last two terms are always singular. The only way that we may obtain bounded 
curvature is to make those two terms cancel between themselves. This requires 

r=\ (52) 

which is exactly the correct scaling for a correction coming from shear viscosity. 
Moreover, we have to fine tune the coefficient A to Ii36i 



A = 2^23"? (53) 
which is the value corresponding to the value of the shear viscosity to entropy ratio 

s Ak 

In this way we reconfirmed, in a fully nonlinear setting ||36l, the value of shear 
viscosity computed at fixed temperature f5|. 

The above analysis can be repeated for the NNLO correction ||38]| with the final 
result for £(t): 



, , 1 2 1 l+21og2 1 
2z3? ^ 12v3 



' Various steps of this calculation were done in 1351 . 1361 and 1371 . 
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The coefficient of the NNLO term involves 2" order transport coefficients. It is at 
this stage that the pathologies of Fefferman-Graham coordinates surface, leading 
to a leftover logarithmic singularity in the scaling limit of the curvature (appear- 
ing at NNNLO in the metric, which is necessary for obtaining the coefficients of 
e(T) at one order lower). The singularity was found to be persistent and not asso- 
ciated with truncating other fields of ten dimensional supergravity |18|. Its origin 
has been explained in detail in |19| (see also |20, 21 1) and can be associated with 
the singular transformation between coordinates regular at the horizon (Eddington- 
Finkelstein) and the Fefferman-Graham ones, coupled with performing an expan- 
sion of the geometry w.r.t. those coordinates. In order to proceed further, which is 
however rather impractical analytically, one would have to perform the analysis in 
Eddington-Finkelstein coordinates (a result in this direction, the NNNLO term in 
e(T) is given in 131]). 



9 Interlude: hydrodynamics redux 

In the above, we have adopted a very agnostic attitude towards the expected dy- 
namics ruling the time evolution of the energy-momentum tensor of the J/ = 4 
plasma system. We did not assume that one could parametrize the T^v in terms of 
such quantities as flow velocity and energy/pressure. We started off from the most 
general Tjjv consistent with the assumed symmetries. By proceeding in this way we 
have an option of describing dynamics which does not fit at all into a hydrodynam- 
ical language. This is in fact the main reason for presenting such an approach here, 
as in the remaining part of the lectures we would like to address the dynamics of 
boost-invariant plasma at small proper times where we do not expect hydrodynamic 
description to be a good starting point. 

On the other hand such flexibility has also significant drawbacks. The determi- 
nation of the transport coefficients of hydrodynamics presented above followed by 
first deriving from AdS/CFT the explicit form of e(T). Given that e(T), one could 
ascertain that the leading term is a solution of perfect fluid equations of motion, 
and together with the subleading term is a solution of viscous hydrodynamic equa- 
tions with a specific value of the shear viscosity (we leave this as an exercise for the 
reader). 

It is thus very interesting to obtain directly the hydrodynamic equations from 
AdS/CFT without making the passage through explicit solutions. This task was per- 
formed in 1 17 1 and is presented in detail in the lectures of V. Hubeny at this school 
||4J. For completeness, let us just summarize here the main idea. 

The static black hole geometry presented in section|4]is dual to a plasma at rest, 
which can be described by a flow vector m'' = ( 1 , 0, 0, 0) and an energy density E (or 
equivalently temperature T). By performing a boost one can obtain a dual geometry 
to a uniformly moving plasma with 4- velocity uy- . 
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In Eddington-Finkelstein coordinates it is given explicitly as 

ds^ = —2ufidx^dr — il j u^Uydx^dx^ + r^{rinv + u^Uv)dx^dx^ (56) 

where r — °o corresponds to the boundary, r = T /k is the horizon while r = is 
the position of the singularity (r = 1/zef c.f. section |4]|. This is an exact solution 
of Einstein's equations. Now promote T and to slowly varying functions of the 
boundary Minkowski coordinates. The geometry ( [56| ceases to be a solution of Ein- 
stein's equations and has to be corrected by terms proportional to gradients. These 
correction terms can be determined with the integration constants fixed by the re- 
quirement of nonsingularity at the horizon. 

Now from the corrected geometry one can read off the T^y which is explicitly 
expressed (similarly to the metric) in terms of T, and the gradients of . The 
numerical constants coming from nonsingularity become exactly the transport coef- 
ficients. In this way one obtains 

s," 

perfect fluid 

-{KT^) fiog2rr +2r,7 + (2-iog2) ( + + 



second order hydrodynamics 



The energy-momentum conservation d^T^^ = of such a T^^ is by definition the 
hydrodynamic relativistic Navier-Stokes equation. 

From the above construction we see that the appearance of hydrodynamics in the 
AdS/CFT correspondence is now completely understood. For any solution of the 



hydrodynamic equations T{xP), u'^{xP), the formula (56 1 and its correction terms 
give an explicit dual metric valid to the same order of the derivative expansion. A 
similar analysis was performed later in the Fefferman-Graham coordinates |39|. 

One final thing to note is that the starting point in the above construction is the 
boosted black hole, which means that one assumes that approximately one is dealing 
with an energy-momentum tensor of a hydrodynamic type (i.e. parametrizable by a 
flow velocity and energy density). This does not always need to be the case, as we 
shall see shortly, and then one has to return to an ab-initio analysis of Einstein's 
equations of the type presented in section |7] 



10 Plasma dynamics beyond hydrodynamics 

The appearance of hydrodynamic behaviour of strongly coupled plasma as a conse- 
quence of the AdS/CFT correspondence is certainly very interesting and satisfying 
theoretically, however perhaps the most fascinating feature of AdS/CFT is its abil- 
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ity to address the behaviour of a plasma system very far from equilibrium, where in 
QCD we do not even have a well motivated phenomenological model. 

As an example of a configuration which cannot be described, even in any approx- 
imation, by hydrodynamic methods consider the problem of plasma isotropisation, 
extensively studied in various variations within QCD ll40l 1411 l42l . 

Suppose we have a plasma system uniform in space with anisotropic pressures. In 
weakly coupled gauge theory one could consider a gas of gluons with non isotropic 
momentum distributions, like gaussians with different widths for the different mo- 
mentum components. Then one expects that the pressures would isotropise in time. 
The energy-momentum tensor of such a system would have the form 



/e \ 

Op\\{t) 
pj_{t) 
Vo px(f)/ 



(57) 



Note that such a system cannot be described by any form of (even all-order re- 
summed) viscous hydrodynamics, since by symmetry — (1,0,0,0) and thus has 
vanishing derivatives. So all viscous terms vanish, while the leading term is clearly 
of a different form. However nothing stops us from applying the AdS/CFT analysis 
using Einstein's equations to such a system. This has been first proposed in |43|. 
Subsequently, a numerical study of this system was performed in B4l . 

Another interesting problem, which we will discuss in the remaining part of these 
lectures, is the behaviour of the boost-invariant plasma system considered before but 
now at small proper times. Since the hydrodynamic expansion i 

, , 1 2 1 l+21og2 1 
T3 2234 ' 12v3 T3 

is an expansion in inverse powers of T, it completely breaks down as we approach 
T = 0. Here the situation is more complicated than in the case of uniform isotropisa- 



tion (57 I mentioned above as we expect a mixture of non-equilibrium and hydrody- 
namic behaviour. In fact it is exactly the question of the transition to hydrodynamics, 
and what factors are important in setting the scale of this transition, that is very in- 
teresting in the context of heavy-ion collisions. A related more general issue is the 
observation of thermalization (proposed in Pi31 to be related to a formation of a 
black hole in the bulk) and an analysis of the concrete way in which this scenario is 
reahsed. 



11 Dynamics at small proper time 

For the reasons described above, we will have to deal with the full Einstein's equa- 
tions. From the point of view of hydrodynamics treated as a gradient expansion 
these encompass all orders of viscous hydrodynamics together with an infinite set 
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of higher transport coefficients. But apart from these there is additional information 
contained in the Einstein's equations. Taking the case of a planar black hole as an 
example, all order hydrodynamics may be identified, on the linearized level, with 
the lowest quasinormal mode and its exact dependence on spatial momentum. But 
apart from this lowest mode there is an infinite set of higher quasinormal modes 
which decay exponentially (in the AdS/CFT context, see in particular |46 |). And all 
of these become important in a far from equilibrium situation such as the early time 
dynamics of the boost-invariant flow. 

Here we will again go to the boost invariant setting and repeat the analysis of 
section |7] but now concentrating on the small T regime. We follow the analysis 
of ll47l . We will adopt the same ansatz for the metric (|34]) and solve Einstein's 
equations with the boundary conditions (|36]). 

Before going into the details of this construction, let us comment why we are us- 
ing the Fefferman-Graham system of coordinates. In contrast to the late time expan- 
sion, here we will aim at solving the Einstein's equations to an arbitrary accuracy 
— without performing any kind of scaling limit. Therefore any choice of coordi- 
nates works equally well. Moreover the constraint equations for initial data in the 
Fefferman-Graham system of coordinates are particularly transparent. 

Let us first determine the qualitative behaviour of e(T) at small T. We will do it 
in two ways. 



The absence of a scaling variable 

In the large proper time regime, the structure of Einstein's equations naturally led to 
the introduction of a scaling variable, which reduced the problem to solving ordinary 
differential equations and a subsequent expansion in inverse powers of T. Let us now 
analyze the solution of Einstein's equations at small T from this point of view. 
We again start from the exact power series solution 

a(T,z)=-£(T).z4 + |-^-^|.Z^ + {i6(T)2 + lT£'(T)e(T) + 

and substitute the asymptotics 

e(T)-— forr-^O (60) 



In this way we obtain 
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(61) 



where the terms dominating for small T are rendered in bold. This analysis was first 
done by Kovchegov and Taliotis B9l . who deduced that for generic s the dominant 
terms at small T are of the form 



^•/(--^) (62) 



In f49\, the scaling solution was found, but due to its complex branch cut stmcturqj 
Kovchegov and Taliotis argued that the only acceptable solution had s ~ 0, which is 
a very interesting result. 



But if we again look at (61 1, we see that the terms resummed by the scaling 
variable vanish for ,5 = and are no longer dominant. Hence there is no place for a 
scaling variable at small T and for s ~Q one has to reanalyze the Einstein equations 
in order to describe the full solution at T ^ 0. 



The existence of a regular initial condition 

One can reach the same conclusion, as well as some more detailed information, on 
the small T dependence of £(t) assuming that at T = we have a regular initial 
conditiorj^ Recall the expression (59 1 and subsitute T = 0. Firstly, we see that the 



assumption that the metric coefficients are finite leads to a finite limit of £(t) as 
T 0, consistent with 5 = 0. Secondly, the inverse powers of T appearing in the 
higher order terms do not lead to a singularity if and only if e(T) has an expansion 
only in even powers of t: 

e(T) = eo + e2T2 + e4T'^ + --- (63) 

A closer analysis reveals that the coefficients e2„ are uniquely determined, 
through the Einstein's equations, in terms of the coefficients of the initial condition 
for the metric: 

a{T = 0,z)=aQZ^ + a2z''+a4Z* + ... (64) 

In order to complete the analysis of the gravitational setup, we have to analyze what 
are the admissible initial conditions (|64]). Once this is done, one can set up the anal- 



And an additional physical bound on s, see |49| . 
' This is an assumption which may, or may not be realistic for ultraenergetic collisions. We prefer 
to keep the options open and analyze boost invariant flow with regular initial conditions as an 
interesting nonequilibrium dynamical system for its own sake. 
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ysis of the system by evolving the geometry from (|64]i using Einstein's equations 
and read off e(T) from the metric. This can be done either numerically solving Ein- 
stein's equations f4T\, or analytically by expressing the coefficients £2,, directly in 
terms of the coefficients of the initial condition a2n BTI . 



The classification of possible initial conditions 

As is well known, in General Relativity, the initial conditions cannot be arbitrary but 
have to satisfy some nonlinear constraint equations. This causes the gravitational 
initial value problem to be quite nontrivial in general. Fortunately, for the case of 
the T = hypersurface in the Fefferman-Graham coordinates, the constraints can be 
solved exactly. 

Let us denote by E^p, the components of Einstein's equations written in the form 

Eap^Rali+'^gali=0 (65) 

Then the constraints are contained in equations /Sfj^ and E„. Denoting flo(z) = fl(T = 
Q,z) etc. we get at once 

ao{z)=bo{z) do = bo = CQ = (66) 

and the only remaining costraint is the single nonlinear equation 

«o +c'i + \ {a'o f + ^ (^o)' - ^ («o + c'o) - (67) 

Let us note an extremely surprising feature of the above equation. At the linearized 
level, it has a trivial solution qq (z) = — cq (z) . So one may expect that for infinitesimal 
ao, the function cq would be also very small and only slightly deformed from — ao 
by taking into account the effect of the nonlinear terms. It turns out, however, that 
this is never true, and the nonlinearity always causes a blowup of the solution for 
some finite z- To see this introduce v(z^) — ^ao(z) and similarly w(z^) for cq. Then 
the constraint equation ( [67] l takes the simple form 

v' + w' + v^ + = (68) 

Now it is easy to see, that there does not exist an everywhere bounded (v = w = 
at infinity) solutions of the constraint equations. To this end it is enough to integrate 
(|68| to get 

poQ poQ r-ao 

0=/ {v' + w')+ {v^ + w^)= (y2 + w2) (69) 
Jo Jo Jo 

Therefore, at T = with boost invariant symmetry, gravity leads to inherently non- 
linear dynamics — a linearized regime does not exist at all! 
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It is not difficult to impose conditions on v and w so that the blowup is not a 
curvature singularity — at least stays finite there (see ATI for details). Moreover 



one can solve the constraint ( 68 1 analytically. Indeed, defining v+ = —w — v, v_ 
w — V, 

2v\-vl (70) 

solves ( [68] l for any v+. Let us conclude with an example of a simple particular 
solution of the initial value constraints: 

ao {z) —bo{z)—2 log cos az^ cq{z) =2 log cosh az^ (71) 

The huge range of initial conditions that can be imposed is in fact quite natural. 
On the gauge theory side we may also expect to be able to prepare an initial state 
with the same energy density in a multitude of ways. Let us contrast this with the 
large T expansion (58 1 of e(T), which only depends on a single scal^^ In other 



words, once the dominant asymptotics of £(t) is known, all subleading powerlike 
terms are uniquely determined. 

The physical interpretation of this difference is quite clear. We expect dissipative 
effects to wash out differences in initial conditions leaving only a single scale (un- 
der the present symmetry assumptions) governing the large proper time expansion 
of £(t). On the gravitational side, these effects can be understood as nonlinear gen- 
eralizations of higher quasinormal modes which die off exponentially (see ll50l for 
some analysis along these lines in the boost invariant setting). 



An analysis of some aspects of the small proper time behaviour 

ofe{x) 

Once we have the allowed initial conditions at T = 0, we need to solve Einstein's 
equations with these initial data. Then, as explained in section |4] we may read off 
£(t) from the solution of Einstein's equations. Since we do not have a scaling vari- 
able at our disposal we have to do it exactly. The ideal way to proceed would be 
to solve Einstein's equations numerically. This study is currently under way ||48ll . 
The route followed in [47] was to solve these equations for the metric in a power 
series in z and T, obtaining a power series expression for e(T) up to the order t'"*' 
for some initial conditions. A drawback of the above method is that the power series 
in question has a finite radius of convergence, necessitating the use of Pade approx- 
imations as an extrapolation method. Below we will present some analysis of these 
extrapolated profiles ||47| . We have to emphasize, however, that one would require 
a real numerical solution of Einstein's equation to be sure of all the details. 

First let us discuss the transition to hydrodynamics. One possibility of quantify- 
ing this is by considering an 'effective exponent' of the power law dependence of 



In pS} this scale has been set to unity, but may be reinstated unambigously by dimensional 
analysis. 
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£(t) defined as 

-T-floge(T) (72) 

Initially it is zero, and it should rise up to 4/3 for late time expansion. In order 
to evaluate it we plot a Pade approximant (of constant large T asymptotics) of the 
expression (|72]). The result, for the initial condition ( [7T| i is shown in Fig. 2. 




We see that it definitely crosses i = 1 of free streaming and moves upwards. 
However, to be sure that it would reach 4/3 a numerical solution for £(t) would be 
needed. 

Now assuming the late time exponent 4/3, we may perform a Pade approxi- 
mation of e(T) with this asymptotics to see the profiles of e(T) for a set of initial 
conditions. Example plots are shown in Fig. 3. 




Fig. 3 Pade resummed profiles of e{T) for a set of initial conditions 1471 . 

The main interest of the knowledge of the exact profile of e(T) for various initial 
conditions is that then we would be able to study the transition to hydrodynamics 
and its dependence on various features of the initial conditions. Since the Pade ex- 
trapolation introduces some serious systematic uncertainties, we refrain from doing 
so until we will have at our disposal a direct numerical solution of Einstein equations 
for these various initial conditions. Apart from the reasons mentioned above, the nu- 
merical solution would also allow to analyze the nature of the apparent singularity in 
the initial data and, more importantly, analyze the geometry for the presence of ap- 
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parent (dynamical) horizons, relevant for the thermodynamic interpretation. Some 
of these issues are currently under investigation 1481 . 

Before we close this last part of the lectures let us note a complementary nu- 
merical investigation of boost invariant plasma in |51 1. The setup of ifSTl was dif- 
ferent from the one presented in these lectures in that the gauge theory metric was 
perturbed in a boost invariant way at some T '-^ To > and then set again to flat 
Minkowski. The metric perturbation produced a change in the geometry which in- 
duced a boost invariant flow on the boundary. The main results observed in lISTl 
were a transition to hydrodynamics and a formation of an apparent horizon which 
moved in from infinity. 

Another related worl|^ was ll52l . were a perturbation by a boundary scalar source 
induced a black hole formation in the bulk. 

Looking at all these examples, one sees that Einstein's equations, through the 
AdS/CFT correspondence, have the potential of describing a multitude of far from 
equilibrium strongly coupled gauge theory phenomenae. It is however clear that the 
majority of problems remain still unsolved. 



12 Conclusions 

In these lectures we have described an approach using the AdS/CFT correspondence 
as a tool for studying real time dynamics of strongly coupled gauge theory plasma. 
The basic theoretical tool is the possibility of translating, in a completely explicit 
and constructive way, between the spacetime energy-momentum tensor character- 
izing the gauge theory configuration in question and the five-dimensional metric 
of the dual geometry. Then one uses the fact, that at strong coupling, the dynam- 
ics of the dual geometry is governed by Einstein's equations (with a cosmological 
constant following from the full 10 dimensional supergravity solution). A further 
dynamical input is the requirement of the absence of naked singularities in the grav- 
ity background. This restricts very strongly the allowed spacetime profiles of the 
gauge theory energy momentum tensor, and consequently the possible gauge theory 
dynamics. 

Using these methods, one may establish the appearance of nonlinear hydrody- 
namics, as exhibited here in the form of near perfect fluid dynamics in the large 
proper time asymptotics of boost invariant plasma expansion. Moreover, Einstein 
equations together with the nonsingularity criterion unambigously predict viscous 
first- and higher-order deviations from perfect fluid dynamics with specific values of 
the transport coefficients appropriate to the case of the ,y/ — 4 SYM theory studied 
here. 

Let us note, that we arrived at these results without presupposing anything even 
about the general form of the gauge theory energy momentum tensor like the pres- 
ence of something like a flow velocity etc. This flexibility comes from the fact 



' This time not in the boost-invariant setting discussed here. 
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that Einstein's equations on the string side of the AdS/CFT correspondence require 
only that the gauge theory energy-momentum tensor is conserved and traceless 
(throughout these lectures we are deahng exclusively with the conformal case of 
jY = \ SYM). Therefore one can use the same techniques to address the question 
of far from equilibrium dynamics where hydrodynamics cannot be used as a start- 
ing point of approximation. We exhibited an example of such a study by describing 
some aspects of the behaviour of the boost invariant plasma expansion in the region 
of small proper time. 

It should be obvious that the majority of questions conceming far from equilib- 
rium dynamics of strongly coupled plasma remain still unanswered even for the case 
of jV = 4 plasma. Once more details are understood, it would be very interesting 
to address similar problems in more compUcated versions of the AdS/CFT corre- 
spondence involving gauge theories which might be closer to QCD. Apart from 
this direct 'application driven' motivation, the study of such time-dependent sys- 
tems leads to fascinating interrelations with General Relativity. On the one hand, it 
provides a new setting for investigating such GR concepts as dynamical apparent 
horizons, black hole formation, providing these GR phenomenae with new inter- 
pretation. On the other hand, the well developed technology of numerical relativity 
might be applied to learn more about the properties of far from equilibrium strongly 
coupled gauge theory systems. 

Last but not least, let us note that the gravity backgrounds obtained as dual de- 
scriptions of evolving plasma systems can very well by themselves form a scene 
for the AdS/CFT correspondence enabling one to study various kind of physics 
questions. In this way one may study how the expanding plasma system influences 
properties of mesons, Wilson loops, various correlation functions. Of course, due to 
the time-dependent nature of those backgrounds this may be quite difficult to do in 
practice, but the possibility of doing so is certainly very interesting and the results 
may be rewarding. 
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Topics not covered in the main text — a brief guide to the literature 

In this appendix we would like to give pointers to the literature on other develop- 
ments related to the approach presented in these lectures. 

• Leading a' corrections (i.e. beyond strong coupling in J/ = 4 SYM theory) to 
the transport coefficients were computed using the boost invariant flow Il53ll54ll . 
This task involved using a' corrected Einstein's equations. 

• Beyond ,y\^ — 4 SYM. A class of general conformal field theories parametrized 
by higher curvature terms in the dual gravitational action was considered in ||55]| . 

• Extension to hydrodynamics with conserved charge(s) was considered in ||56l 
|57]|58l. Electric-magnetic equilibration at large proper times was found in ll56ll . 
In addition, dilaton driven hydrodynamics was considered in a general way 
in ll59l and the onset of turbulence was observed ll60l . 

• Lower dimensional examples. The case of a 1h-1 dimensional conformal field 
theory allows for an explicit exact dual gravitational solutions ||6T1 |62l . Other 
investigations in different number of dimensions were performed in 1631 . 

• Various (exact) solutions for ,y\^ = 4 gauge theory in curved and possibly time- 
dependent backgrounds were obtained 1641 l65l \66\. The exact solutions which 
were found do not involve viscosity effects. An exact gravitational description af 
a shearless flow in ^ = 4 in flat space was obtained in |67 1. 

• Further properties of solutions with boost invariant symmetries were studied 

in ESllelllTollTllCll. 

• Physics in the expanding plasma. Fundamental flavours were introduced (through 
D7 brane embeddings) into the late proper time boost invariant geometry Ii73l . 
diffusion constant was evaluated ll74ll . drag force was computed |75 1. In addition, 
various physical questions were addressed in the case of the shock wave solutions 
described in section]?] Deep Inelastic Scattering (DIS) was analyzed 1176117711781 
|79J as well as heavy quark energy loss 180. 81J . 
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